APS/123-QED 



Role of soft-iron impellers on the mode selection in the VKS dynamo experiment 
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A crucial point for the understanding of the von-Karman-Sodium (VKS) dynamo experiment is 
the influence of soft-iron impellers. We present numerical simulations of a VKS-like dynamo with 
a localized permeability distribution that resembles the shape of the flow driving impellers. It is 
shown that the presence of soft-iron material essentially determines the dynamo process in the VKS 
experiment. An axisymmetric magnetic field mode can be explained by the combined action of the 
soft-iron disk and a rather small a-effect parametrizing the induction effects of unresolved small 
scale flow fluctuations. 
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Homogenous dynamo action, i.e. the transformation 
of the kinetic energy of a flowing conducting fluid into 
magnetic energy, accounts for the generation of plane- 
tary, stellar and galactic magnetic fields. The last decade 
has seen tremendous progress in the experimental inves- 
tigation of dynamo action [l[. By exploring parameter 
regions that are hardly accessible to numerical simula- 
tions, such experiments are essential for a better under- 
standing of geo- and astrophysical magnetic fields. So far 
the threshold of magnetic field self-excitation has been 
crossed only in three liquid sodium experiments in Riga, 
Karlsruhe, and Cadarache. In the latter one, the von- 
Karman-Sodium (VKS) experiment, a turbulent flow of 
liquid sodium is driven by two counter-rotating impellers 
located at the opposing end caps of a cylindrical domain, 
and field generation sets in when the magnetic Reynolds 
number exceeds the critical value of Rm c w 32 [2|. Be- 
sides of the rich dynamical behavior, including periodic 
and irregular field reversals that occur when the impellers 
rotate with different frequencies |3[ , the structure of the 
measured field attracts attention because of its high de- 
gree of axisymmetry jij. This geometry is in contrast 
to all numerical predictions based on axisymmetric flows 
from which a magnetic field with an azimuthal wavenum- 
ber m = 1 (B tx cost/?) was expected 5-7]. Furthermore, 
dynamo action is only obtained when the impellers are 
made of soft-iron with a relative permeability of the or- 
der of /i r ~ 100. So far, there is neither a satisfying 
explanation of the observed axisymmetric field, nor of 
the failure of the dynamo in case that non-magnetic im- 
pellers are being used. This means, in turn, that the 
very working principle of the VKS dynamo is not under- 
stood at all. A first attempt to explain the axisymmet- 
ric mode in the VKS experiment was recently made in 
terms of a turbulent aw-model 8. 9]. The consideration 
of turbulence seems imperative since the kinetic Reynolds 
number is of the order Re ~10 6 , and flow and field ex- 
hibit strong fluctuations that are of the same order as 
the mean values. The model allows to circumvent the 
restrictions imposed by Cowling's theorem which forbids 



axisymmetric dynamo action (Io| . Indeed, it was shown 
in |11| that an axisymmetric VKS-like flow supported by 
an a-cffcct is able to excite an axisymmetric field. How- 
ever, the magnitude of the a-effect which is necessary 
to explain the dominant axial dipole turned out to be 
several times larger than the rms value of the turbulent 
velocity. Therefore, this simple aw-model seems not ca- 
pable of explaining the VKS dynamo mechanism. But 
what about the role of the soft-iron impellers? Origi- 
nally, their use was motivated by the hope to shield off 
the field in the bulk of the cylinder from the lid layers 
behind the impellers. These layers, and particularily any 
rotating flow therein, had been shown to be extremely 
detrimental for dynamo action [7]. First attempts to 
investigate the role of the high permeability impellers 
utilized simplifying boundary conditions that enforce a 
vanishing tangential field at the top and the bottom of 
the cylindrical domain which mimics the field behaviour 
in case of end caps with infinite /x r 0, [l2| . A similar ap- 
proach was followed in [l3| with focus on the impact of 
an additional non- axisymmetric velocity component with 
azimuthal wave number (m = 8) that describes a hypo- 
thetic vortex structure between adjacent impeller blades. 
Comparable to the aw-model with a localized a-effect 
flit ' an axisymmetric field was obtained when the am- 
plitude of the distortion is of the same order as the mean 
Sow. This means, however, that a similar problem as for 
the aw-model arises, because the assumed existence of 
highly correlated, large amplitude non-axisymmetric mo- 
tions could not be confirmed in water experiments [lij ]. 
In this Letter, we will treat the rotating soft-iron im- 
pellers in a more realistic manner than just to assume 
vanishing tangential field conditions at the end caps. We 
consider dynamo action caused by an axisymmetric mean 
flow u which might be supported by an a-effect. The 
large scale magnetic flux density B is governed by the 
mean field induction equation 

dB „ / „ 1 „ B 
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^ocr(r) /i r (r). 
Here a denotes the conductivity, fj,o the permeability of 
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vacuum and /x r the relative permeability. An a-effect 
might be caused either by a gradient of the turbulence 
intensity or by helical fluid motions predominantly in- 
between the eight blades of the impeller disks. In the 
latter case, the (m = 8) flow modulation is considered 
as a small scale variation whose averaged influence is 
parametrized by an axisymmetric a-source term 0, fll| . 
From Eq. ([I]) it follows that a high (even if only localized) 
permeability /x r reduces the effective magnetic diffusivity 
rj cB = (afMopf)- 1 , where [if = V~ x J ^(r)d 3 r and V 
denotes the cylinder volume. In addition, the terms on 
the r.h.s. of Eq. (Q} that involve gradients of a or /i r 
are able to couple toroidal and poloidal field components 
in a way that possibly allows exponentially growing so- 
lutions for B. A paradigm for this type of dynamo has 
been presented by Busse and Wicht [15| who showed that 
even a constant flow over an infinite plate with spatially 
varying conductivity is able to sustain a dynamo. Inter- 
estingly, the structure of the field is then determined by 
a constant part modulated by a contribution with the 
wave number of the conductivity variation. A similar 
effect results from the inhomogeneous permeability dis- 
tribution in the VKS experiment, and it is worthwhile 
to check if this is already sufficient to explain a signifi- 
cant axisymmetric contribution. However, a purely ax- 
isymmetric field cannot be expected as such a solution 
is forbidden even in case of space- and time dependent 
permeability /conductivity distributions (16|. 

Eq. (ED) is solved numerically applying a finite volume 
method [17J . Note that at the experimental temperature 
of ~ I20°C, the differences of the electrical conductivity 
of iron and liquid sodium are not very important so that 
we focus on a non-uniform permeability distribution. In 
the grid based approach the implementation of spatially 
varying material coefficients is straightforward and only 
requires the application of a simple averaging directive 
for the reconstruction of the local (discretized) value of 
/i r that is seen by a field component at a certain grid cell. 
Although the scheme is easy to implement, the simula- 
tions are still demanding because large gradients of /i r 
must be resolved. To keep the computational time in a 
reasonable limit, the presented results are restricted to 
a maximum value of /Lt r = 60. This is in rough agree- 
ment with recent measurements of the permeability of a 
VKS-like impeller that gave values of fi T ~ 65 ... 70 [l8[ . 

A suitable mean flow that describes the velocity distri- 
bution in the VKS-experiment is represented by the so 
called MND-flow [II], 

u r (r,z) = -Tr/Hcos(2%z/H)r(l -r) 2 (l + 2r), 
Utp(r,z) = 4er(l - r) sin (ttz/H) , (2) 
u z (r,z) = (l-r)(l + r-5r 2 )sm(2nz/H), 

where H = 1.8 denotes the distance between both im- 
peller disks and e describes the relation between toroidal 
and poloidal component of the velocity. Here, e = 0.7259 



which is the optimum value for obtaining dynamo action, 
at least when the mean flow alone is considered. The en- 
tire model setup is visualized in the left panel of Fig. [Q 
The velocity field ^ is only applied in the region between 
the two impellers. The flow active region is restricted to 
a cylinder of radius R = I (corresponding to 20.5 cm in 
the experiment) which is embedded in a layer of stagnant 
fluid with a thickness of OAR. This so-called side layer 
significantly reduces Rm c [3] . A lid layer is added behind 
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FIG. 1: Left panel: Structure of the prescribed axisymmet- 
ric velocity field. The color coded pattern represents the az- 
imuthal velocity and the arrows show the poloidal velocity 
field. The black solid lines represent the shape of the impeller 
system (disk and blades). Right panel: assumed permeability 
distribution; isolevel at fj, T = 35 where the peak value is given 
by fi r = 40. In the fluid region fi r is equal to unity. 

each impeller disk. Within these lid layers a purely ro- 
tating flow is assumed, modeled by a linear interpolation 
of the azimuthal velocity at the outer side of the impeller 
disk towards to zero at the end cap of the cylindrical do- 
main. The presence of lid- and side layers reduces the 
influence of the magnetic boundary conditions so that, 
for the sake of simplicity, vanishing tangential field com- 
ponents are imposed at the outer surface of the compu- 
tational domain. The solid black lines in the left panel 
of Fig. [Q represent isolevels of the permeability distribu- 
tion that mimics the impellers with radius r — R and 
height h = 0.2. The impeller system is composed of two 
fractions: the disk, which is modeled by an axisymmet- 
ric scalar field with constant /i r > 1, and an azimuthally 
varying contribution located adjacent to the inner side 
of each disk. This part represents the system of eight 
blades and is modeled as an azimuthal rotating scalar 
field oc cos(8</5 + u>t) in a way that the minimum value 
(located between the blades) corresponds to fi T = 1 (the 
permeability of the fluid) and the peak value corresponds 
to the disk permeability. Here, u> denotes the angular ve- 
locity of the impeller that is necessary to drive a flow 
with a given Rm (defined as Rm = Rfio<j\U ma _ x \) . The 
3d structure of the resulting permeability distribution is 
shown in the right panel of Fig. [T] Note that the impellers 
are modeled only by the permeability distribution and no 



3 



particular flow boundary conditions are enforced on the 
(assumed) interface between impeller and fluid. 

Initially, simulations have been performed without any 
a-effect. After a transient phase that lasts approximately 
one diffusion time r ~ R 2 /r) cS , the system settles in an 
eigenstate characterized by exponential growth, respec- 
tively decay, \B\ ~ e 7 * where 7 denotes the field am- 
plitude growth rate. Fig. [5] shows the growth rates for 
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FIG. 2: Growth rates of the field amplitude (without a-effect) 
Solid curves represent the results for the axisymmetric mode 
(m — 0) and dashed curves show the (m = l)-mode. The 
inset plot presents the growth rates for /j, t — 1 and fi r = 60 in 
dependence on Rm. 

the axisymmetric field mode (to = 0, solid curves) and 
the simplest non- axisymmetric mode B occost/? (to = 1, 
dashed curves). Even at moderate values of fi r the 
(to = 0)-growth rate is considerably shifted towards the 
dynamo threshold (but remains still negative). The re- 
duction of this decay rate is roughly in agreement with 
the reduction of r) e g . For large permeabilities the growth 
rate of the axisymmetric mode saturates slightly below 
the dynamo threshold. In the experimental relevant in- 
terval (RmSs50) the (to = 0) growth rate dominates and 
remains nearly independent of the applied Rm. Contrary 
to the axisymmetric mode, the (to = l)-mode is hardly 
affected by the permeability distribution (for Rm <; 15), 
and for sufficiently large Rm, the (to = l)-mode domi- 
nates again, resulting in dynamo action around Rm c 
60. Preliminary simulations applying only an axisym- 
metric distribution for the permeability (modeling the 
impeller disc without blades) exhibit a reduction of the 
(to = l)-mode to Rm c rs 50 whereas the axisymmetric 
field mode is hardly influenced (for Rm ^ 50). 

It is important to note that so far no growing ax- 
isymmetric solutions were obtained. One might have ex- 
pected a sort of Bussc-Wicht dynamo with an axisym- 
metric component modulated by a (unobserved) higher 
azimuthal wavenumber due to the azimuthally varying 
blade system. But evidently, this mechanism alone is 



not sufficient to explain the occurrence of the dominant 
axisymmetric contribution in the VKS experiment. In 
the following we will see that an axisymmetric growing 
solution can be obtained if an additional small a-effect 
is invoked. For the sake of simplicity we will assume 
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FIG. 3: Growth rates of the axisymmetric mode for Rm = 30 
in dependence of a prescribed homogenous a-effect (ot vv ). In 
all cases the (m=l)-mode lies below the growth rates of the 
axisymmetric solution. Additionally, the hypothetical case 
with u po \=u r e r + u z e z =0 has been studied for fi r = 20. The 
inset plot shows the scaling behavior of the critical value for 
a. vv in dependence on fj^ s . 

a uniform distribution of a vv . This component is re- 
sponsible for the generation of the poloidal field from 
the toroidal field whereas the generation of the toroidal 
field essentially is determined by shear so that the other 
components, a rr and a ZZl are assumed to be less im- 
portant. Because the high permeability impellers have 
already shifted the growth rate of the (to = 0)-mode 
very close to zero, it is not surprising that a small a- 
effect is now sufficient to excite a growing axisymmetric 
field even at moderate values of the impeller permeabil- 
ity (see Fig. [3]). For larger /i r the dependence of the 
growth rates on the parameter a is rather delicate and 
even a small variation in a results in large changes of 
the corresponding growth rates. Growing axisymmetric 
fields are generated for both signs of a, although they 
are obtained a bit more easily for positive values of a. 
The critical value of a vv that is necessary for growing 
axisymmetric fields is strongly decreasing for increasing 
/i r (see inset plot in Fig. [3]). In physical units the crit- 
ical value for /j, r = 40 corresponds to |a c | ~ 0.15 ms _1 
which is much lower than the critical value reported in 
[111 ] where |a crlt | « 16 ms _1 was obtained (for /i r = 1 and 
Rm = 32). The local maxima for /i r = 20 (blue curve) 
and fx v = 40 (black curve) close to a = correspond 
to stationary solutions whereas a transition to oscillat- 
ing solutions occurs approaching the adjacent minima 
which might be the cause for the local suppression of 
the growth rate. Interestingly, the hypothetical omission 
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of the poloidal velocity component (u po \=u r e r +u z e z =0) 
works even in favor of the dynamo (see red curve in 
Fig. [3]). This might explain the intriguing experimen- 
tal observation that, despite a careful optimization of 
the impeller blades for one rotation direction, a change 
of this blade orientation does not alter significantly the 
critical Rm As long as the a-effect close to the im- 
pellers is strong enough, poloidal field generation is suf- 
ficient to sustain the dynamo cycle. The resulting mag- 
netic field is strongly concentrated in the impeller region 
(see Fig. For larger Rm, the field geometry is de- 

termined by the well known (m = 1) dominated field 
structure (Fig. 4b) and dynamo action predominantly 
takes place in the bulk of the cylinder which explains the 
impassivity of the (m = l)-mode towards the high /i r 
disks located close to the end caps. Note that in case 




(a) (b) (c) 



FIG. 4: Distribution of magnetic energy density (blue col- 
ored isosurfaces at 10% of the maximum value) and geometric 
structure of the magnetic field (red streamlines) for fi T = 20: 
(a) Rm = 30, a — —1.5, growing axisymmetric mode; (b) 
Rm = 70, a — 0, growing (m = l)-mode; (c) Rm = 30, a — 0, 
decaying axisymmetric mode. Note the (m = 8) modulation 
induced by the impeller blades. 

of a vanishing a-effect the (decaying) solution is domi- 
nated by an axisymmetric toroidal field and the poloidal 
field generation is not sufficient to close the dynamo cy- 
cle (Fig. Hh). At this point, a deeper investigation re- 
quires a more detailed modeling of the a-effect on the 
basis of water experiments. Direct numerical simulations 
of forced turbulence are far away from realistic param- 
eters so that reliable models for the a-distribution can 
hardly be extracted from DNS models. Preliminary sim- 
ulations with a different spatial distribution of a suggest 
that the extreme facilitation of the axisymmetric mode 
by the soft-iron disks is a generic effect which takes place 
as long as sufficient poloidal field generation (due to a) 
occurs nearby the impellers. 

The dynamo process is strongly influenced by the re- 
duction of rf s . However, the decrease of rf s is not suffi- 
cient to explain the results because simulations with the 
permeability variation replaced by a conductivity varia- 
tion of the same magnitude (and space/time dependence) 
exhibit a significant deviation in the magnitude of the 
growth rates (although the tendency remains similar) . 

Our results clearly demonstrate that the common sim- 



plified treatment of the high permeability impellers, by 
demanding vanishing tangential field components at the 
interface, is not justified. In the relevant range of Rm, 
the high permeability works selectively in favor of the 
(rn = 0)-mode while the (m = l)-mode is only slightly 
affected. To initiate a growing (m = 0)-mode, some a- 
effect has still to be invoked, although the scaling behav- 
ior with fi r indicates that its magnitude is very small. 
The homogenous a-distribution should be regarded as 
the simplest example how the soft iron disks facilitate 
growing axisymmetric solutions at reasonable parameter 
values. A complementary and more detailed approach 
will have to consider a non-axisymmetric flow variation 
in terms of azimuthally drifting equatorial vortices that 
have been observed in water experiments reported by 
20]. Similar as in the a-model with high permeability 
blades this should lead to a mixed mode solution of the 
magnetic field where a higher azimuthal mode is present. 
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